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FOREWORD > 


The advent of the nuclear submarine with its unlimited 
cruising range has increased the urgency of acquiring and an- 
alyzing data required to assess the ocean environment. This 
report presents one example of the type of analysis required 
to understand some of the characteristics of random-type data. 
Such data are exemplified by continuous recordings of sub~ 
marine motions and relative flow across the deck of the 
hovering or slowly moving submarine, and digital recordings 

of the height of the sea surface above the submarine. The 

fact that so many interrelated variables can be sensed and 
recorded simultaneously is an indication of the submarine's 
potential as an oceanographic research platform. In addition, 
such simultaneous recording procedures make the data amenable 
to even more sophisticated analysis; namely, the estimation of 
cross-spectral densities connecting two continuously distributed 


random variables. 


Bo Go Sill? 
Rear Admiral, U. S. Navy 
Hydrographer 
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1, INTRODUCTION 


The submarine, USS REDFIN (SS-272), was utilized as a working platform by the 
Submarine Systems Section, U.S. Navy Hydrographic Office, to obtain digital re- 
cordings of surface wave motion and analog recordings of ship and fluid motions. This 


report is concerned with the spectral analysis of these recordings made while the 
REDFIN was hovering. 


One of the objectives of the Submarine Systems Section is to collect and analyze 
data needed to study environmental effects on submarine motion, and the characteristics 
of the relative water motion in the proximity of the submarine. Such data, recorded 
while the submarine was hovering or underway at the lowest possible speed, have been 
analyzed by computing power spectra of the recordings. 


Data presented in this report include composite graphs of spectra of surface wave 
motion (computed by digital estimation formulae), spectra of the submarine's roll angle, 
and spectra of the transverse and logitudinal relative flow velocity of water across the 
submarine's deck (obtained by an analog computer). Examples of spectra of in situ 
pressure fluctuations, estimated by both digital and analog methods, are presented. 
Differences in approach embodied in the digital versus analog method of power spectra 
estimation are discussed, and specific examples are presented to illustrate points of 
technique required for comparable resolution and statistical stability. Finally, spectra 
data presented in this report are summarized by presenting mean values of motion in 
tabular form. 


The objective of this particular test was to study environmental effects on the 
motion of the submarine. During this test period, wave conditions were generally mild 
with the result that some data were recorded near the sensitivities of the instruments. 
Thus, energy spectra and related statistical parameters are representative of conditions 
encountered under lower sea states. 


This report includes most of the data recorded on REDFIN Cruisel! during the 
period 28 October to 5 November 1959. During this time, the REDFIN was hovering 
at a keel depth of 100 feet, and operating in an area centered about 180 miles east of 
Wilmington, North Carolina. 


Hydrographic Office personnel participating in this test were Messrs. Q. H. 
Carlson (Senior Oceanographer), A. G. Alexiou, and D. E. Tidrick. 


Il, INSTRUMENTATION AND RECORDING PROCEDURES 


A. Instrumentation 


Wave, ship, and fluid motion observations discussed in this report were made from 
aboard the REDFIN. Instrumentation utilized and variables measured are included in 
the listing that follows. Locations of these instruments on the submarine are shown in 
Figure 1}. 


Pitch and Roll Angle David Taylor Model Basin stable 
platform and potentiometer 

Relative Flow Velocity Litton Electromagnetic Flow Meter 

In Situ Pressure Fluctuations Wiancko Pressure Measuring System 
(U.S. Navy Mine Defense Laboratory) 

Ship's Speed Ship's Electromagnetic Log 

Ship's Course Ship's Gyrocompass 

Depth (pressure) Vibrotron and Bourns potentiometer 


in the pressure gauge line 


Surface Wave Profile Westinghouse Sonic Surface Scanner 


Figure 2 illustrates a portion of a typical recording taken aboard the REDFIN. 
Most data were recorded on both magnetic and paper tape. Surface wave and in situ 
pressure fluctuation data, however, were recorded on magnetic tape only. A com- 
plete description of the instrumentation and measurement and recording procedures is 
given in the Hydrographic Office Technical Report No. 91 (Reference 10). 


Nearly all pitch angle recordings, heave acceleration, and sway acceleration 
were practically flat and will not be discussed further in this report. Thus, the only 
submarine motion of Interest is roll angle. Some cross-flow data with very low 
amplitudes were deleted. 


Although attempts were made to provide analog recordings of in situ pressure 
fluctuations measured by the Wiancko pressure gauge, most of the collected data were 
deleted because the amplitude of the recording's envelope frequently exceeded the 
dynamic range of the pressure gauge. When this occurred, the relief valve in the gauge 
would open, resulting in intermittent jumps in the continuous recordings. Attempts to 
filter these pressure data before spectral analysis were not too successful, and most of 
the spectra of the Wiancko recordings were not considered valid. 
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Surface wave data were recorded by the Sonic Surface Scanner in the form of 
pulses on magnetic tape. Since this method may prove to be a useful way of obtaining ~ 
wave data in the open ocean while submarines are on routine patrol, the theory and 
method of this technique are reviewed briefly in Appendix A. 


B. Recording Procedure 


The recording procedure was as follows: The REDFIN submerged to 100 feet (keel 
depth) and attempted to hold this depth under zero ground speed. Recordings of one- 
half hour duration were made under these conditions at various headings and wave 
conditions. Test conditions are summarized in Table 1. Exact headings were difficult 
to maintain; headings, relative to wave directions, shown in Table 1 are the average 
values during each run. Each recording was given two numbers for identification. For 
example, Run 3-1 is the first recording taken on Run 3. Except for surface wave data 
and in situ pressure fluctuations, each run was recorded on paper tape for visual anal- 
ysis and monitoring of data recorded at sea, and, also, on magnetic tape for electronic 
analysis. 


Most surface wave data in Table 1 were obtained by "hindcasting" from synoptic 
weather maps (Reference 9). There were a few visual observations, but these observa~ 
tions are considered unreliable, limited in many cases by darkness. The general diffi- 
culties of making visual wave observations at sea are discussed in Reference 11. Hind- 
cast wave data provide an independent source of wave height, period, and direction. 
These are particularly useful as a check on the accuracy of visual observations made 
aboard a submarine, in the event no surface data are available from the Sonic Surface 
Scanner . 


Hi. THEORY AND METHOD OF POWER SPECTRUM ANALYSIS 


A. Spectrum Analysis 


A preliminary analysis of most of the variables discussed in this report was pre- 
sented in Reference 8. This report contains mean and maximum values of the peak-to- 
peak roll angles, and peak-to-peak cross-flow velocities. However, motion charac- 
teristics of such variables (including surface wave motion) can be more completely 
described by computing energy spectra of the recordings. Analog and digital recordings 
of motion data described in this report are assumed to be samples from quasi-stationary 
Gaussian processes whose mean values are zero. Such processes are completely described 
by their energy (power) spectra. 


Suppose part of a given recording of the process occurs as shown in Figure 3. 
Although such functions are random, amplitudes of the process (determined from a line 
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connecting the zeros of the process) can be given statistically by a probability density 
function. If the dominant part of the random process occurs over a narrow Interval of 
frequencies, amplitudes shown in Figure 3 can be given by the Rayleigh probability 
distribution: 


2 
Bell! 
E 


F (x) dx = 


e€ 


dx (1) 


Thus, the probability that a particular amplitude a, will fall between x and x + dx 
is given by equation 1. The first moment about the origin of equation 1 gives the re- 
lationship between the mean wave amplitude a and the quantity E: 


<= 0-866 VE- (2) 


The quantity E is given by the second moment about the origin of the Rayleigh 


distribution: 
eo 
2x3 ee 
Be eae crm maak (3) 
6 


E_.9 
Bare (4) 


where o” is the variance of the process. Now E is related to the energy spectrum of 
the random process by: 


= 2 2 
= LAY dx and — = LAG (5) 


° 
where (ROY is called the energy density. From the energy (power) spectra and 
related statistical parameters, the process can be completely described for practical 
purposes. The term "power spectra" is freely used by aero- and hydro-dynamicists 
among others, although much of the work dealing with the energy spectra of random 
processes originated in the field of communications engineering (Reference 5). 


10 


Recordings of submarine motion (roll angle) and fluid motions (transverse and longi- 
tudinal relative flow velocity) illustrated in Figure 2, may be visualized by considering 
them frequency band-limited random signals continually oscillating above and below an 
equilibrium level of zero signal. Very often (especially in the case of wave and ship 
motions), the spectrum of such motions is characterized by a well-defined peak value 
which is skewed toward lower values of frequency, The area under the spectrum curve 
is related to mean amplitude values of the motion, and from the theoretical work of 
Longuett—Higgins (Reference 2) et al, the following relations connecting mean wave 
amplitude a to the means of the one-third highest (a,) and the one-tenth highest wave 
amplitudes (a, ) were found: = 


il 3 

10 
an i 
10 = 1.800 me (6) 
a (o] 


These amplitude relationships were verified on actual wave height data by J. K. A. 
Watters (Reference 12), who found the following mean wave height relationships H. 
from a sample of wave records: 


1 
& 

= 1.94 a2 = | ,58 (7) 
A 


From these we conclude that wave height values are given approximately by doubling 
the amplitude values. Finally, since: a@ = 0.866VE then 


ia =v =) 177 ee (8) 
from which we get: 
H, = 2.33 /E | (9) 
3 
H, = 3,60 JE (10) 
0 


VW 


Thus, equations 8, 9 and 10 hold for mean values of roll angle and relative flow 
velocity, since recordings of roll angle and flow velocity are assumed to be samples 
from the same type stationary Gaussian random process. 


B. Digital Estimation of Power Spectra 


The practical problem of estimating energy spectra has received widespread 
attention and can be viewed from two distinct but practically equivalent points of 
view. One method of estimating the energy spectrum of a stationary Gaussian random 
process is to sample the continuous (analog) recording of the process at discrete inter- 
vals of time, obtain an estimate of the autocorrelation function, and then obtain raw 
or unsmoothed estimates of the energy spectrum by Fourier transformation of the auto- 
correlation function. The power spectrum Is then obtained by suitably smoothing the 
raw or uncorrected spectral estimates. Since the study of such processes is the study 
of the sample variance of the process, ft is essential that confidence intervals be given 
for each spectrum. The theory and procedures of digital estimation of power spectra 
are given in References 1, 3, and 6. 


If the process can be "quantized" or digitized at the transducer or in the recording 
process, much of the manual data reduction can be eliminated. For example, the Sonic 
Surface Scanner measures instantaneous height of the sea surface above the transducers 
located on the hovering submarine's deck, and these data are recorded directly in 
digital form on magnetic tape. These values are played back, put on punched cards, 
and programmed through a high speed computer to obtain the estimated spectra of the 
surface wave heights. The method of recording and playback is outlined briefly in 
Appendix A, 


As an example, consider Figure 4, which shows the power spectrum obtained from a 
one-half hour recording from transducer No. 1, located on the submarine's bow (Fig. 1). 
The dashed lines indicate the 90-percent confidence intervals (See Section C). For 
comparison, the spectrum of a simultaneous one-half hour recording from transducer No. 
8 (located near the stern) is shown in Figure 5 with its 90-percent confidence intervals. 
Figure 6 shows the good agreement between the two measured spectra. 


Amplitudes of submarine motion during this recording were low. The heave accel- 
eration recording, for example, was practically flat. However, at higher sea states and 
lower wave frequencies, a submarine will move appreciably in response to surface wave 
motion. For example, heave displacement would be in phase with a pure swell, whose 
relative heading was 90-degrees, with the result that the fraction of the true recorded 
wave amplitude decreases with decreasing frequency (See Appendix A for a discussion of 
expected measurement errors). In general, surface wave spectra estimated from record- 
ings of the Sonic Surface Scanner, will have to be corrected for the hovering submarine's 
motion. Furthermore, if the submarine is underway, an additional correction will be 
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necessary for its speed and relative heading. 


Results of power spectra of surface wave motion obtained on the global circum- 
navigational cruise of the USS TRITON (SSRN-586)* indicate that such spectra can be 
measured successfully at speeds of 6 knots. At 6 knots, the TRITON's speed was equal 
to the phase speed of a 2-second wave. Thus, spectra given in terms of the frequency 
of encounter could be transformed in terms of the frequency relative to fixed coordi- 
nates. 


C. Analog Estimation of Power Spectra 


The digital method of estimating power spectra may not be practical If a large 
number of recordings are to be analyzed. This being the case, the recommended pro- 
cedure would be to record analog data on a multichannel tape recorder and compute 
energy spectra by commercially available wave analyzers. This was the method used 
to analyze ship- and fluid-motion data presented in this report (References 4 and 7). 
Most of the analog spectra of submarine and fluid motions presented in this report were 
computed at the David Taylor Model Basin. The Hydrographic Office has recently 
acquired the electronic components required to estimate power spectra shown in the 
block diagram of Figure 7. An example of an analog spectrum of submarine roll angle 
estimated by the system at the Hydrographic Office is shown in Figure 8. 


The analysis procedure is as follows: The random process to be analyzed is recorded 
aboard the submarine on magnetic tape at a tape speed of 1 7/8 inches per second. Each 
one-half hour run recorded at this speed is played on a 14-channel Ampex playback unit 
(Fig. 7) at 60 inches per second while being re-recorded on a 14-channel Precision 
recorder. The random signals are re-recorded also on an 8-channel Sanborn recorder so 
that the quality of the recordings can be inspected, and an indication of expected am- 
plitudes of the random signals can be obtained. This results in a compression of 1/32 
the original recorded tape length; the frequency components in the random signal are 
all increased by a factor of 32. The compressed recording is formed into a continuous 
loop and threaded into a special loop attachment. The Precision recorder is utilized 
again to drive the loop attachment at a speed of 30 inches per second for an additional 
frequency increase of 16 times and a total frequency increase of 32 x 16 = 512 times 
the original frequencies. Thus, a frequency of 0.10 cycles per second appears on the 
analyzer as a frequency of 51.2 cycles per second. The reason for analyzing in this 
frequency range is that beat frequency analyzers commonly used can easily accommodate 
frequencies at higher ranges, but could not handle gravity wave frequencies (e.g. 0.03 
to 0.30 cycles per second). 


“U.S. HYDROGRAPHIC OFFICE, Oceanographic data report for the global circum- 
navigational cruise of the USS TRITON (SSRN-586). 16 February to 10 May 1960, by 
NR. Mabry (in press). 
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The random signal output is fed into the analyzer system consisting of an oscillator, 
analyzer, and power integrator. These three units transform the random time signal into 
a continuous plot of the average energy density versus frequency, which is traced out on 
a Moseley X-Y Recorder. This system also provides for estimation of the integrated spec- 
tra. The block diagram of Figure 7 illustrates an integrated spectrum curve superimposed 
on the energy density curve. The small inset curves in the lower right hand corner of 
the X-Y Recorder represent the spectral outputs of the filter used to resolve the oridinate 
scales of the graphs of average energy density and the integrated spectrum. Details and 
circuit diagrams are given in "Instruction Booklets for TP-625 Wave Analyzer System, " 
Technical Products Company, Los Angeles, California. A discussion of the problems 
associated with filter bandwidth, loop periods, time constants, and frequency scanning 
rates is given in Reference 4. 


Figure 8 is an example of the power spectrum of a submarine roll angle recording 
and the power spectrum of the filter used in the analysis. This illustrates a typical 
power spectrum of the filter used to calibrate the ordinate scale of the energy density 
graphs. The broken line curve is the digital spectrum of the same recording which is 
superimposed for comparison, The digital spectrum superimposed is one of four spectra 
shown on Figure 9 (See Section D on comparison of spectral estimates) . 


This analog spectrum was estimated by the wave analyzer in the Hydrographic 
Office, and all the digital spectra contained in this report were estimated by a 
Burroughs 205 Electronic Computer, also located in the Hydrographic Office. 


Such spectral estimates are distributed with a chi square distribution with 2N degrees 
of freedom, where N is the number of elemental bands covered by the filter. If the 
spectral density is not a fast changing function of frequency interval equal to the band- 
width of the filter, then N is approximately equal to the effective bandwidth of thé 
filter divided by the elemental bandwidth. The elemental bandwidth is 1/T, where T is 
the loop period, and the effective bandwidth used was 6 cycles per second (Reference 4). 
Thus, for the analog spectra in this report: 


effective bandwidth 
elemental bandwidth 


Number of degrees of freedom = 2N = 2 (11) 


Confidence intervals for the estimated spectra can be determined from a table of the 
chi square distribution. 


The 90-percent confidence intervals for the digital spectra are given by computing 
the number of degrees of freedom according to: 


2 
Number of degrees of freedom = TR (12) 
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where Tn is the length of the record in seconds, m is the number of lags, and At is the 
sampling interval. For example, the spectral estimates in Figures 4 and 5 are distrib- 
cuted with a chi square distribution with 70 degrees of freedom. The 90-percent confi- 
dence intervals can be determined from a table of chi square distribution. For 70 
degrees of freedom, the true spectral density Pe (f) is bounded by the estimated spectral 
density Pe (f) according to: 


PC) Pe (f) 


a 30 < PS (f) < wOu76 or 0.77P_ (FKP (F)<1.32P, (F) (13) 


Thus, if a single observed estimate with 70 degrees of freedom is observed to be 


5 ft 2/sec™!, then we have 7 Uae confidence that ne true Igng-run value lies 
between 5/1. 30 = 3.85 ft 2/sec™! and 5/0.76 = 6.58 ft 2heec 


D. Some Considerations in Estimating Spectra 
1. Resolution and Statistical Stability 


In estimating power spectra, whether by digital or by analog methods, a 
choice must be made between a high resolution of the spectral estimates over a chosen 
frequency band and the statistical stability of the spectral estimate in that band. The 
two considerations are mutually opposed in that fulfillment of one implies some sacrifice 
of the other. For example: 


a. Width of the "spectral window" determines the resolution of any peaks in 
the spectrum. Now, if the spectral window (i.e., the filter in the wave analyzer) is 
chosen too narrow, too much detail is obtained, and a truly significant peak may be 
overlooked in the resultant "blurring." However, if the spectral window is made too 
wide, the resulting spectral estimate may be too smooth, and here again some physically 
significant hump in the spectrum may be smoothed over. 


b. The statistical stability of the estimate depends on the width of the spectral 
window and the length of the sample, When high resolution is desired, fewer Fourier 
components are averaged over the narrow spectral window, and the spectral estimate will 
have wider fluctuations over the totality of ensemble estimates. Conversely, less reso- 
lution with a wider filter gives a value which is expected to deviate less from the true 
but unknown spectrum. 


2. Aliasing 
a. If asignal g(t) contains no frequencies above f cycles per second, it is 


completely determined by giving its ordinate by a series of points spaced Sr seconds 
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apart. For example, if waves with frequencies greater than 0.5 cycles per second are 
absent, the wave record can be sampled every second. However, if appreciable wave 
energy is present at frequencies above f cycles per second, this energy will contribute 
(by aliasing) to the spectrum of the lower frequency waves. Thus, the sampling inter- 
val must be determined by the high frequency cut-off. 


b. The range of wave frequencies shown on the abscissas of Figures 4, 5, 10, 
and 11 represent the most important part of the spectral energy density associated with 
surface wave heights. Similarly, the range of frequencies in Figure 9 are associated 
with the dominant rolling motion of the submarine. Now, if a random process is char- 
acterized by oscillations about a slowly varying centerline, the spectral estimates at 
frequencies near zero may have appreciable values. However, since this spectral 
energy is not associated with frequencies in the gravity wave spectrum, which are the 
only frequencies of Interest in this report, the graph is arbitrarily cut off at the lowest 
important frequency. Similarly, spectral densities at the high frequencies were not 
plotted since spectral densities decreased rapidly as the Nyquist frequency was 
approached, 


3. Comparison of Spectral Estimates 


When the estimated analog spectrum of a process is to be compared to a 
digital estimate of the same spectrum, the equivalent spectral windows must be the 
same. Figure 9 illustrates 4 digital estimates of the analog spectrum shown In Figure 
8, The analog spectrum in Figure 8 had 58 degrees of freedom, and the four digital 
spectra had 31, 42, 63, and 125 degrees of freedom. Therefore, the digital spectra 
with 63 degrees of freedom were superimposed on Figure 8 for comparison. 


IV, POWER SPECTRA OF SURFACE WAVE, SUBMARINE, 
AND FLUID MOTIONS 


A. General 


Wave, ship, and fluid motion data analyzed in this report were taken while the 
REDFIN was executing a series of maneuvers which included hovering at different 
relative headings. During each hovering run, the REDFIN attempted to hold a keel 
depth of 100 feet at a fixed relative heading while making recordings of surface wave 
height (Sonic Scanner), roll angle, transverse relative flow velocity, longitudinal 
relative flow velocity, and in situ pressure fluctuations. Each recording was at least 
one-half hour in duration. 
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Table 1 is a summary of test conditions for Runs 3, 4, 5, and 6 of Cruise Il. The 
first recording began with Run 3. Wave conditions were generally mild except for the 
first four recordings of Run 3. On four of the runs (Runs 4-1, 5-1, 5-4, and 6-3), large 
variations in depth occurred because of temporary loss of buoyancy control. However, 
such large depth variations were of very short duration, and the submarine generally 
was able to return to its original hovering position In about 7 minutes. All runs where 
this occurred were continued for 30 minutes after depth was regained, and only this 
portion of the record was analyzed. 


Submarine motions and relative flow velocity of water across the submarine's deck 
occur as a result of surface wave motion. Since motion characteristics of submarines 
will be required under various sea states, It is essential that some estimate of the surface 
wave spectrum be available. Such spectra can be obtained by hindcasting from the 
synoptic sea-level surface weather maps by methods described in H. O. Pub No. 603 
(Reference 9). However, the exact functional form of the Pierson-Neumann spectrum 
and its dependence on such parameters as wind speed, wind duration, and fetch length 
are presently being Improved. Hence, It is essential to obtain a series of surface wave 
spectra over the open ocean so that differences between measured and predicted spectra 
can be reconciled. 


Although the Sonic Surface Scanner method js still in the experimental stage, it is 
a promising method of obtaining series of surface wave records in the open ocean. Sub- 
ject to corrections for the motion of the submarine, it should be possible to obtain reli- 
able spectral estimates of surface wave heights from such records. 


B. Surface Wave Height Spectra - 


Surfeice wave height spectra for Runs 4 and 6 are shown in Figures 10 and 11, . 
respectively. The sonic scanner was inoperative during Runs 3 and 5. Figures 10 and 
11 indicate that wave conditions were stationary during each of these series of runs. 
No corrections have been applied to these surface wave spectra, since the amount of 
wave energy associated with high period waves is rather low. Table 2 presents the 
90-percent confidence intervals for these spectra. 


Each of the submarine and fluid motion spectra shown in the composite graphs 
(Figs. 12 to 21) was obtained in the same form as depicted in Figure 8. For ease of 
comparison and tdentification, however, a smooth curve was reconstructed from each 
analog spectrum. 


C. Roll Angle Spectra 


Figures 12 through 15 present analog energy spectra of roll angles for Runs 3 to 6. 
Each figure contains 4 spectra corresponding to a different relative heading. Each roll 
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angle spectrum has a typical peak, occurring between 16 and 20 seconds, which is 
associated with the natural rolling period of the REDFIN. The natural rolling period 
of the REDFIN is about 18 seconds. Secondary peaks in Run 3, occurring at periods 
lower than 10 seconds, are associated with the lower period waves observed during 
this run (Table 1). Table 3 presents the 90-percent confidence intervals for all analog 
spectra discussed in this report. 


TABLE 2. NINETY PERCENT CONFIDENCE INTERVALS FOR 
DIGITAL SPECTRA OF SURFACE WAVE HEIGHTS 


Degrees of 

Run No. Freedom Lower Limit Upper Limit 
Sa ee eS see i Sy ln he lal eres 

4-] 80 .78 Pe (f) 1.27 PE (f) 

4-2 80 .78 (NA7-7A 

4-3 80 .78 27. 

4.4 80 .78 Te27 

6-1 60 74 1.39 

6-2 60 74 1.39 

6-3 66 .76 Wee 7/ 

6-4 68 76 lmao 


D. Transverse Relative Flow Velocity Spectra 


Figures 16 and 17 present energy spectra of transverse relative flow velocity for 
Runs 4 and 6, respectively. Figure 16 presents 5 spectra while Figure 17 presents 2 
spectra, each corresponding to different relative headings. Spectra for Runs 3 and 5 
were deleted because the amount of energy associated with frequencies below the 
lower cut-off value (.03 cycles per second) of the wave analyzer made the determina- 
tion of parameters extremely doubtful. In fact, the excessive energy produced some 
spectra which appeared as "white" noise (flat spectrum) from 0.03 to 0.20 cycles per 
second (5- to 33-second periods). This problem will require further study, 


The two higher peaks in Figure 16, occurring in the first two recordings of Run 4, 
ore associated with periods of 8.3 to 9.0 seconds, corresponding to dominant periods 
of surface wave motion. Peaks at higher period values are associated with the cross- 
flow energy induced by the rolling motion of the submarine, which tends to roll at its 
resonant period. Similar peaks occur in Figure 17, with the peaks centered at lower 
periods associated with the prevailing lower period wave motion. 
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TABLE 3. NINETY PERCENT CONFIDENCE INTERVALS 
FOR ANALOG SPECTRA 


Degrees of 
Run No. Freedom Lower Limit Upper Limit 
3-1 46 we Pe(f) 1.47 Pe(f) 
3-2 35 68. 1.61 
3-3 31 .68 1.67 
3-4 32 .67 U ole? 
4-] 58 5H) 1.39 
4-2 49 72 1,45 
4-3 4] 71 1.49 
4-4 48 52 1.45 
A-5 38 ./0 1.54 
5-1 29 .68 1.64 
5-2 Sy) 74 1.37 
5-3 37 ./0 ieoZ 
5-4 25 65 eae 
6-1 28 .68 1.67 
6-2 AO 71 1.49 
6-3 37 ./0 1.52 
6-4 24 66 2 


E. Longitudinal Relative Flow Velocity Spectra 


Figures 18 through 21 present energy spectra of longitudinal relative flow velocity 
for Runs 3 to 6, respectively. Each figure, except Figure 20, contains four spectra 
corresponding to a different relative heading. 


Relative headings for the first two recordings in Figure 18 are designated with 
respect to sea because these relative headings associated with sea waves had values 
that could result in larger longitudinal flow components, even though swell was greater 
in height and period bandwidth (Table 1). In the last two recordings, relative headings 
with respect to swell are considered most important, although the choice is arbitrary. 
The extreme peak in Run 3-1, however, probably was caused by the joint contribution 
of both sea and swell. 


In Run 3-4, both sea and swell components were from abeam of the submarine, 


and this results in the lowest spectral values throughout the run because the longitudinal 
component is smallest. 
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FIGURE 13, ROLL ANGLE SPECTRA FOR DIFFERENT RELATIVE HEADINGS UNDER HOVERING CONDITIONS, RUN 4 
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Similarly, in Runs 4 and 5 (Figs. 19 and 20), spectral energy values Increased 
during those runs due to the favorable relative headings. During Run 6 (Fig. 21), 
peaks definitely shifted to lower period values, reflecting the low period wave values 
at the surface. 


V. SUMMARY OF RESULTS 


Table 4 compares values derived from the measured surface wave height spectrum 
to similar values obtained by the hindcasting methods outlined in H. O, Pub. No. 
603. Significant wave heights obtained by hindcasting are lower than values obtained 
from the measured digital wave spectra. In addition, the highest periods from hind- 
casting are 7.8 to 8.8 seconds while the highest periods from the measured spectra are 
12.0 to 14.4 seconds. Wind speeds are plotted on the Northern Hemisphere synoptic 
sea-level maps in increments of 5 nautical miles per hour. For fully-developed condi- 
tions, a significant wave height of 3.2 feet and a wave period band of 1.5 to 7.8 
seconds occur when the wind speed is 14 knots; at 18 knots (fully developed) a signifi- 
cant wave height of 6.1 feet and wave perlod band of 2.5 to 10.0 seconds would occur. 
Thus, a 4-knot increase in wind speed (an increase of about 5 miles per hour) results in 
a significant wave height increase of 3 feet and an increase of 2.2 seconds at the high 
end of the wave-period band. In addition, high period swell of very low amplitude 
from many different sources, whose contribution could be detected in a measured wave 
spectrum, could never be hindcasted from the present oceanic weather maps. Thus, 
data presented in Table 4, from two different sources, appear to be fairly consistent. 


TABLE 4. SUMMARY OF RESULTS DERIVED FROM SURFACE WAVE SPECTRA 
(DIGITAL) AND COMPARISON TO HINDCAST WAVE DATA 


Values from the Measured Values from Hindcast 


Nave Spe C 


H 1/3 T -Band 
(Ft) (Sec) 


4-2 5.9 2.5-14.4 3 1,5-7.8 
4-3 5.4 2.5-14.4 3 1,5-7.8 
4-4 5.6 2.5-14.4 3 1,5-7.8 
4-5 Dink 2.5-14.4 4 2.0-8.3 
6-1] 6.3 2.4-12.0 5 2.0-8.8 
6-2 5.6 2.4-12.0 5 2.0-8.8 
6-3 6.0 2.4-12.0 5 2.0-8 .8 
6-4 6.7 2.4-12.0 5 2.0-8 .8 
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Table 5 presents a summary of some of the Important parameters derived from the 
analysis of roll angle spectra. The first column gives the run number of each recording; 
the second, an indication of the success (or failure) of the submarine in attempting to 
hold to a fixed keel depth of 100 feet; and the next three, mean values of the peak-to- 
peak roll angles in accordance with equations 8, 9, and 10. The last three columns 
give mean values of peak-to-peak roll angles computed manually for comparison. The 
agreement in most cases is good. 


Tables 6 and 7 summarize relative flow velocity data. Columns 3, 4, and 5 are 
the mean values derived from spectrum relationships in accordance with equations 8, 
9, and 10. Columns 6, 7, and 8 are the mean values obtained manually. The agree- 
ment between the spectrum derived mean values and those computed manually is good 
in most cases. 


Vi. CONCLUSIONS 


This report Is concerned with the measurement and analysis of certain variables 
applicable to oceanographic problems; the chief analysis tool is the energy spectrum. 
Specific conclusions are: 


1. The Sonic Surface Scanner method appears to be a reliable method of 
obtaining surface wave records in the open ocean. 


2. Parameters derived from the estimated spectrum obtained by the Sonic Surface 
Scanner compare reasonably well to similar parameters obtained by hindcasting. 


3. As expected, the largest amplitudes of roll angle occur near the REDFIN's 
resonant frequency. A secondary maximum, associated with frequencies of maximum 
spectral energy of the surface wave motion, is observed in almost every roll angle 
spectrum. 


4, Spectra of transverse relative flow velocities tend to exhibit two distinct 
peaks: one at frequencies near the resonant frequency of the REDFIN, and a secondary 
peak is well correlated with the surface wave motion. 


5. Spectra of longitudinal relative flow velocities exhibit maxima whose frequen- 


cies are well-correlated with the frequency of maximum spectral energy of the surface 
wave motion. 
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6, Although spectra of roll angle and cross~flow velocities were presented with 
relative heading as parameter, no consistent relationship is Indicated by these spectra 
data. This is attributed largely to the low sea state values which persisted throughout 
the test except on Run 3. 


Vil. FUTURE TEST AND ANALYSIS PLANS 


1. Data discussed in this report are representative of lower sea states. To 
correctly assess the influence of environmental conditions on submarine motions, 
additional data are required during State 4 or 5 sea conditions. Such data should 
be taken under conditions of head seas (O-degrees relative heading), abeam seas 
(90-degrees), and following seas (180-degrees), in order to study the characteristics 
of motion at several different relative headings. 


2. This report presents some examples of "auto" spectra, or the resulting 
transformation when the wave or ship motion record is autocorrelated (correlated 
with itself). The extension of the concept leads to cross-spectra, or the transforma- 
tion resulting when simultaneous recordings of two different random processes are 
cross-correlated (correlated with each other). A random function of time representing 
a stationary Gaussian process (for example, sea surface wave motion) could be con- 
sidered the input or forcing function. The output or response function would be sub- 
marine motion. Now if frequencies exist in the input such that an excitation response 
analogous to "resonance" is possible in the output, then such a frequency interdepend- 
ence is clearly displayed by cross-spectrum analysis. Thus, the use of cross-spectrum 
analysis provides a useful tool in the study of cause and effect relationships between 
two random processes. 


3. Recent communications between leading oceanographers indicate that the 
array of nine sonic scanners may be used to obtain the directional spectrum of the 
natural seaway under a variety of sea conditions. Such spectra would be estimated 
by cross-spectrum analysis. Sea surface wave motion is a two-dimensional phenomenon, 
and energy of wind-generated surface waves has been studied as a function of frequency 
alone (the notable exception being Stereo Wave Observation Project). 
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APPENDIX A = SURFACE WAVE HEIGHT DATA OBTAINED BY THE 
SONIC SURFACE SCANNER 


A. Digital Recording and Playback of Sonic Surface Scanner Data 


An example of a sea surface profile obtained by Sonic Surface Scanner Is presented 
in Reference 8. The Sonic Surface Scanner consists of ten transducers mounted on the 
deck, These are spaced 33 feet apart (Fig. 1). The transducers are mounted face 
upward and are essentlally inverted echo sounders. Each transducer In turn measures 
distance between Itself and the sea surface directly above in a 3° cone. By sequencing 
from one transducer to the next aft along the deck of the submarine, a wave profile, 
comparable to the length of the submarine, ts obtained. The complete sequence takes 
approximately 0.7 seconds. The Sonic Scanner was designed to provide an aes) 
to the nearest foot. 


The output of all ten transducers are available, but the primary interest In this 
report is the power spectra of the output of one transducer. A recording from one 
transducer is, therefore, the instantaneous height of the sea surface above the submarine 
as a function of time. This is given in terms of the two-way travel time required for a 
sound wave front to travel from the transducer to the sea surface and return. 


Assume the velocity of sound in ocean waters to be 5000 ft/sec (the resulting error 
in this assumption is very small). Travel time per foot can be expressed: 


sec ./ft. = 0.2 x 1073 sec ./ft. = 0.2 millisec./ft. 


1 
5000 


Since we are concerned with a two-way travel time (the time required for the 
sound beam to go from the submarine transducer to the sea surface and return), the 
effect is the same as If the travel time were doubled. Thus, the effective travel time 
per foot is: 


= 0.4 mlillisec ./ft. 
For convenience In counting, t is expressed as: 
= 10 millisec./25 ft. 
In this way, changes tn the height of sea surface above the submarine are measured 
in terms of the two-way travel time. These changing values of sea surface height above 


the submarine are recorded as pulses on an AM channel of a tape recorder every 0.7 
seconds. 
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Assume the sequence of values recorded at Transducer No. 1 Is desired. Appro- 
priate pulse signals are monitored on an osciHoscope (Fig. Al), and the time Interval 
between outgoing pulse and return echo is counted on the time Interval counter. 
Time values are then printed out on a digital recorder. As a check on the quality of 
the digital printout, an Esterline-Angus Recording milliameter is used to record the 
analog output from the digital recorder, which appears as a staircase function and is 
smoothed out by the slow response time of the Esterline-Angus Recording milliameter . 


Normally, the counter displays time In milliseconds, and this Is proportional to 
the sea height above the submarine. For convenience in computation, it Is desirable 
to have the print-out In feet directly. This Is accomplished by introducing a 25 ke 
signal to the external time base of the time Interval counter. The output of the 
printer is then given directly in feet and tenths of feet. These values are coded sub-~ 
sequently on punched cards for high speed computation of power spectra on the 
Datatron. 


B. Effect of Submarine Motion on Surface Wave Measurements 


Sonic Surface Scanner wave measurements at a keel depth of 100 feet and made 
under lower sea state conditions will not be affected seriously by the submarine 
moving in response to the subsurface wave motion. However, as the sea surface wave 
height grows under the action of the wind, or if low frequency swell propagates over 
the hovering submarine, the low frequency (high period) components will increase 
the amplitudes of motion of the submarine. In effect, wave measurements will be made 
with respect to a moving reference system. 


Now, if water depth is great, the attenuation factor of a simple harmonic wave in 
the neighborhood of the hovering submarine fs given by: 


, 2 
alum -2nZ. agi: 
KS Ruker for L = On 

where: L = wave length 


T = wave period 
zZ depth to the submarine 


Ht} 


Depth in the operating area off Wilmington, North Carolina, was 250 fathoms, 
and we will consider KS for Z = 75 feet and 100 feet. These two depths represent 


average depth to the Sonic Surface Scanner transducers and the average keel depth 
for Sonic Scanner spectra data shown In'Figures 10 and 11. Expressing K, asa | 
Cine ~ | ee H Z 
function of frequency f for fixed depth Z: 
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Figure A2 presents the wave height attenuation as a function of frequency alone 
for depths of 75 and 100 feet. 


Thus, for a low frequency swell of .05 cycles per second (20-second period) only 
20 percent of the surface wave height witl be attenuated at 75 feet, and only 26 per- 
cent will be attenuated at 100 feet. If the submarine behaves as a particle, the 
heave amplitude of the neutrally buoyant, hovering submarine should be something 
like a simple harmonic motion In phase with the simple harmonic wave at the surface 
(neglecting for the moment, all other components of submarine motion). In any event, 
a correction must be applied to either the orlginal data or possibly to the spectrum. 


In general, the submarine will execute motion in all six degrees of freedom. The 
effect of pitching and rolling will tend to Increase the path length over which the sound 
beam has fo travel; this tends to make the recorded wave motion appear higher than the 
actual wave motion, which overcomes some of the difficulty mentioned In the preceding 
paragraph. 


Both these difficulties can be ameliorated !f submarine data are recorded at a depth 
of 200 or even 300 feet. However, with Increased depth, there will be some "skipping" 
of data since, occasionally, a return echo will be too weak to gate the time Interval 
counter. 
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APPENDIX B ~- POWER SPECTRAL ESTIMATES 
OF IN SITU PRESSURE FLUCTUATIONS 


A. General 


One of the varlables the Submarine Systems Section attempted to record was 
pressure fluctuation, as measured by a Wiancko differential pressure gauge. This 
gauge was mounted on the sail with the diaphragm parallel to the deck of the sub- 
marine (Fig. 1). Thus, when the submarine hovered at a keel depth of 100 feet, the 
gauge's rubber diaphragm was actually 60 feet below the surface. 


Total pressure change occurring at the diaphragm of the differential pressure 
gauge of the hovering submarine will be composed of, (1) pressure change due to 
surface wave motion, (2) linear motions of the submarine (heave, sway, and surge), 
and (3) angular motions (roll, pitch, and yaw). 


Pressure changes at the gauge due to surface wave motion alone are independent 
of the submarine's relative heading. Total motion of the submarine, hence total 
pressure change assoclated with the submarine's motion, is definitely a function of 
relative heading. Even under stationary wave conditions, power spectra of the 
Wiancko recordings will, in general, be different for different relative headings. 


B. Spectral Estimates of Filtered Signals 


The problem of recording pressure fluctuations on the hovering REDFIN was 
complicated by the limited dynamic range of the Wiancko differential pressure gauge. 
To prevent damage to the gauge, a relief valve eliminates differential pressure in the 
interior chamber of the gauge when the dynamic range of the gauge Is exceeded. A 
typical in situ pressure recording is shown in Figure B1. The small oscillations are 
pressure changes due to surface wave motion and linear and angular components of the 
submarine motion. The very long "period" oscillations are associated with the slowly 
changing dive angle of the hovering submarine. The amplitude ratio of the large to 
small oscillations is about 20 to 1. This highly compressed recording represents approxi - 
mately 8 minutes of recording. Many of the Wiancko pressure recordings were charac~ 
terized by breaks or "jumps" which occurred when the amplitude of the long-period 
oscillations became too large. When this happened, the relief valve would relieve 
the differential pressure momentarily, and the recording would begin anew at the 
center of the chart. 


Power spectra of the recording shown in Figure B1 were estimated by the wave 
analyzer system at the Hydrographic Office and are shown in Figure B2. In Figure B2, 
three different estimates are superimposed to show the effect of filtering. Curve A is 


the spectrum of the unfiltered signal while Curves B and C are the resultant spectra of 
first passing the signal through a Krohn-Hite Band=Pass Filter, Model 320A. Cut-off 
limits were arbitrarily set at 0.08 to 0.5 cycles per second (periods of 2 to 12.5 
seconds) for Curve B and 0.09 to 0.5 cycles per second (periods of 2 to 11.1 seconds) 
for Curve C, Filter response curves for the frequency band used to obtain Curve B 
are shown in Figure B3. It should be noted that the shape of the filter response curve 
allows a fraction of the energy associated with frequencies less than 0.08 cycles per 
second to pass, but blocks off the very low frequencies. 
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FIGURE Bl, COMPRESSED RE-RECORDING OF IN SITU PRESSURE FLUCTUATIONS 


A digital spectral estimate of this same run was obtained by re-recording the 
filtered run with pass-band limits set at 0.05 to 0.5 cycles per second. In Figure B4, 
the dashed curve is the power spectrum of this filtered recording; the solid curve, 
surface power spectrum; and the dashed-dot curve, the residual power spectrum at 60 
feet obtained by attenuating the surface power spectrum according to: 


2 
a Zo feet (- = LAY, 


where Z and L are given in Appendix A. 


2 


urface 


The difference between the in situ power spectrum and the residual spectrum at 
60 feet due to surface wave motion alone is attributed to pressure fluctuations associ- 
ated with the submarine motion; i.e. roll angle and cross-flow. 
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